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Support  for  Advanced  Imaging  of  Premixed  Turbulent  Combustion  Processes 


1)  Statement  of  Problem  Studied 

For  a  broad  range  of  turbulent  flow  conditions  the  effect  of  flow  turbulence  on  premixed 
flames  is  to  wrinkle  and  distort  the  flame  sheet  thereby  increasing  combustion  rates  per 
unit  volume  -  the  combustion  intensity.  In  previous  ARO  sponsored  research  we 
developed  crossed-plane  imaging  methods  for  the  study  of  reaction  sheet  wrinkling 
(crossed-plane  laser  tomography)  and  changes  in  reaction  sheet  thermal  structure 
(crossed-plane  Rayleigh  imaging)  due  to  turbulence.  In  the  present  research  we  combine 
crossed-plane  laser  tomography  (CPLT)  with  stereo  particle  image  velocimetry  (SPIV)  to 
measure  flamelet  wrinkling,  flamelet  speeds  and  the  reactant  velocity  field  in  order  to 
relate  flamelet  properties  to  the  turbulent  velocity  field  in  the  reactant  flow. 

And  Summary  of  Important  Results 

Using  combined  CPLT  and  SPIV  we  have  obtained  a  large  data  set  for  1 1  different 
flames  stabilized  in  a  V-flame  burner  and  5  different  flames  stabilized  in  a  swirl-flame 
burner.  We  are  in  the  process  of  analyzing  these  data  and  have  concentrated  so  far  on  the 
V-flame  burner  data.  The  most  significant  finding  to  date  is  the  presence  of  negative 
displacement  speeds  and  a  correlation  between  displacement  speed  and  flamelet 
curvature.  Negative  displacement  speeds  have  be  predicted  but  this  is  the  first  time  they 
have  been  observed  in  turbulent  flames.  Details  of  the  measurements  and  results  to  date 
are  given  below. 
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Gouldin  and  S.  S.  Sattler,  AIAA  2004-0976,  presented  at  the  42nd  Aerospace 
Sciences  Meeting,  January  2004,  Reno,  NV. 

(d)  Downstream  Development  of  Premixed  Turbulent  V-Flames,  R.  Beresky,  S.  S. 
Sattler  and  F.  C.  Gouldin,  AIAA  2005-0153,  to  be  presented  at  the  43nd 
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(c)  Support  for  Advanced  Imaging  of  Premixed  Turbulent  Combustion  Processes, 
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6)  Scientific  Progress  and  Accomplishments 

ADVANCED  IMAGING  OF  PREMIXED  TURBULENT  COMBUSTION 
PROCESSES 

F.  C.  Gouldin 

Mechanical  and  Aerospace  Engineering 
Cornell  University,  Ithaca,  NY  14853 

Introduction 

Over  the  last  several  years  at  Cornell  we  have  developed  and  applied  advanced 
imaging  methods  -  crossed-plane  laser  tomography  (CPLT)  [1  -  5],  crossed-plane  laser 
Rayleigh  imaging  [6]  and  stereo  particle  image  velocimetry  (SPIV)  —  for  high-resolution 
studies  of  scalar  and  velocity  fields  in  premixed  turbulent  combustion.  Pulsed  laser 
tomography  with  micron  sized,  silicone  oil  droplets  can  be  used  to  image  the 
intersections  of  650  K  isothermal  surfaces  and  the  laser  illumination  plane.  We  have 
used  CPLT  to  measure  the  instantaneous  orientation  of  these  surfaces  and  to  determine 
the  probability  density  function  of  the  surface  normal  vector,  a  statistical  measure  of 
surface  wrinkling  that  is  used  to  calculate  the  mean  surface  to  volume  ratio  or  the  surface 
density  of  the  650  K  surface.  More  recently,  we  have  combined  CPLT  and  SPIV  [7]  in 
order  to  simultaneously  measure  the  surface  normal,  the  reactant  velocity  field  and  the 
reaction  sheet  (flamelet)  displacement  speed  in  the  laboratory  reference  frame  and 
relative  to  the  local  reactant  gas  velocity.  In  addition  numerical  experiments  have  been 
used  to  evaluate  the  potential  for  using  CPLT  to  measure  the  principle  curvatures  of  the 
wrinkled  reaction  sheet  in  addition  to  the  sheet  normal  vector. 

The  goal  of  our  current  research  is  to  exploit  combined  CPLT  and  SPIV  in  a 
systematic  study  of  premixed  turbulent  combustion  in  a  V-flame  burner  and  a  swirl 
burner.  The  expected  result  of  this  study  will  be  a  unique  database  that  quantifies  both 
reaction  sheet  wrinkling  by  turbulence  (reaction  sheet  normal  vectors,  principle 
curvatures  and  surface  density  data)  and  velocities  (velocity  field  of  reactants, 
displacement  speeds  and  velocity  correlations  at  the  reaction  sheet  surface). 

To  this  end  we  have  made  measurements  of  flamelet  displacement  speed  in  a  methane 
-  air  V-flame,  looked  at  measurements  of  curvature  in  V-flames,  built  a  swirl  burner 
based  on  a  design  by  Robert  Cheng  at  Lawrence  Berkeley  National  Laboratory,  and  taken 
combined  CPLT  and  SPIV  data  on  a  total  of  eleven  turbulent,  premixed  methane-air,  V- 
flames  and  on  five  swirl  burner  flames.  These  developments  are  outlined  in  this  report. 

The  people  working  on  the  research  covered  by  this  report  are  the  PI,  Sandra  Sattler, 
a  Ph.D.  student,  Soo  Ho  Lee,  a  M.  Eng.  student  and  three  undergraduate  students  -  Ian 
Bell,  Ryan  Beresky  and  Joonil  Kwak.  Research  into  curvature  determination  was 
reported  at  the  last  AIAA  Aerospace  Sciences  Meeting.  In  addition  a  manuscript 
reporting  on  previous  work  support  by  ARO  was  submitted  to  Combustion  and  Flame  [6] 
and  has  been  accepted. 
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Method  and  Quantities  Measured 

Combined  Crossed-Plane  Laser  Tomography  and  Stereo  Particle  Imaging  Velocimetry : 
SPIV  allows  one  to  measure  a  2-D  distribution  of  the  3-D  velocity  vector  field.  For 
SPIV,  two  CCD  cameras,  mounted  in  a  Scheimpflug  configuration,  are  focused  from 
different  angles  onto  the  same  plane.  Each  camera  takes  two  exposures  of  laser  light 
scattering  from  seed  particles  in  the  flow  at  two  closely  spaced  points  in  time.  The 
measured  displacements  of  the  seed  particles  between  exposures  are  used  to  determine 
particle  and  thereby  gas  velocities.  A  single  camera  can  be  used  to  obtain  the  in-plane 
velocity  components;  images  from  a  pair  of  cameras  allow  one  to  estimate  the  out-of¬ 
plane  velocity  component  as  well  as  the  in-plane  components.  A  double-pulsed,  dual¬ 
head,  frequency  doubled,  130  mJ/pulse  Nd:YAG  laser  fired  sequentially  in  time  is  the 
light  source  for  these  measurements.  The  beams  from  the  two  laser  heads  are  aligned 
collinearly  and  formed  into  overlapping  sheets  in  the  image  plane  of  the  cameras.  The 
laser  and  image  acquisition  cameras  are  controlled  by  a  data  acquisition  and  analysis 
computer.  A  commercial  SPIV  system  including  laser,  optics,  cameras  and  computer  is 
used  for  the  measurements.  The  camera  and  laser  sheet  configuration  used  for  combined 
SPIV  and  CPLT  is  shown  in  Fig.  1. 

In  CPLT,  images  of  seed  particle  scattering  are  taken  near  simultaneously  from  two 
intersecting,  orthogonal  laser  sheets.  The  seed  particles  are  micron  sized  silicone  oil 
drops  that  evaporate  on  or  near  the  650  K  isothermal  surface  [8].  In  each  image  a  curve 
marking  the  disappearance  of  the  droplets  and  therefore  the  650  K  isotherm  is  identified. 
Tangent  vectors  to  these  curves  projected  onto  the  illumination  planes  are  tangent  vectors 
to  the  650  K  surface.  Where  this  surface  crosses  the  line  of  intersection  between  the  two 
illumination  planes,  two  tangent  vectors  can  be  determined  and  their  cross  product  is 
normal  to  the  650  K.  surface  at  the  instant  the  images  are  taken.  Repeated  laser  firing  and 
image  pair  acquisition  makes  it  possible  to  generate  PDF’s  of  surface  normal  vectors. 

Using  oil  seed  with  combined  SPIV  and  CPLT  one  can  obtain  the  distribution  of 
reactant  velocities  in  a  vertical  plane  and  the  650  K.  surface  normal  vectors  over  a  100 
microsecond  time  period.  These  velocities  and  surface  normal  vectors  are  found  at  the 
three  points  where  the  650  K  surface  crosses  the  three  horizontal  lines  of  intersection 
between  the  laser  illumination  planes.  These  lines  are  marked  a,  b  and  c  in  Fig.  2.  In 
addition,  by  staggering  the  pulse-timings  of  the  tomography  laser  and  the  SPIV  lasers  by 
approximately  1.5  milliseconds,  displacement  speeds  of  the  650  K  surface  in  the 
laboratory  frame  and  relative  to  the  gas  flow  are  determined. 

The  data  obtained  from  combined  CPLT  -  SPIV  measurements  can  be  used  in  several 
different  ways.  Normal  vector  data  can  be  used  to  find  the  650  K  sheet  surface  density, 
an  important  modeling  term  that  is  proportional  to  the  mean  rate  of  product  formation. 
The  PDF  of  the  sheet  normal  vector  is  one  measure  of  the  wrinkling  of  the  sheet  by  the 
turbulence  as  is  the  PDF  of  the  principle  radii  of  curvature.  Reaction  sheet  structure  is 
perturbed  by  the  turbulent  flow  field  through  imposed  flow  strain  and  curvature.  The 
reactant  field  velocity  data  obtained  with  SPIV  and  CPLT  can  be  used  to  estimate 
elements  of  imposed  strain.  Also  of  importance  is  the  displacement  speed  on  the  cold 
side  of  the  reaction  sheet,  which  is  the  difference  between  the  component  of  the  gas 
velocity  on  the  reactant  side  normal  to  the  reaction  sheet  and  the  sheet  speed  in  the 
laboratory  frame. 
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Figure  1.  Schematic  diagram  of  combined  SPIV  and  CPLT  system  showing  the  overall 
layout  of  the  system  including  lasers,  cameras,  optics  and  burner.  The  laser  sheets 
propagate  out  of  the  plane  of  the  paper.  To  avoid  over-exposing  the  cameras  to  laser 
pulses,  the  CPLT  images  are  taken  first.  The  SPIV  images  are  taken  approximately  1.5 
milliseconds  later.  Surface  normal  vectors  are  measured  where  the  650  K  surface  crosses 
the  three  lines  of  intersection,  denoted  by  a,  b,  and  c  in  figure  2,  between  the  three 
illumination  planes  defined  by  laser  sheets  1  and  2  and  the  SPIV  laser  sheet.  Velocity 
vectors  are  obtained  over  the  vertical  SPIV  laser  sheet.  The  field  of  view  of  each  of  the 
four  cameras  is  approximately  35mm  x  35mm.  The  tomography  cameras  have  512x512 
pixels,  while  the  PIV  cameras  have  1360  x  1030  pixels. 


Figure  2.  Schematic  of  combined  CPLT-SPIV  laser-sheet  orientations.  1  and  2  are  the 
perpendicular  CPLT  laser  sheets  and  lie  at  45°  with  respect  to  3,  the  vertical  SPIV  laser 
sheet.  The  sheets  intersect  forming  three  measurement  lines:  a,  b,  and  c.  Flamelet 
normals  and  displacement  speeds  are  measured  along  lines  b  and  c;  flamelet  normals  are 
measured  along  a. 


6 


Results  to  Date 

Data  Taken :  Data  sets  consisting  of  CPLT  and  SPIV  images  have  been  collected  for  1 1 
different  V-flames.  The  flame  parameters  are  listed  in  Table  1,  and  the  location  of  each 
flame  on  the  Peters  regime  diagram  [9]  is  shown  in  Figure  3.  Each  data  realization 
consists  of  two  tomography  images  and  four  SPIV  images;  for  each  flame  condition 
images  from  approximately  2000  realizations  have  been  obtained  for  data  analysis. 
Similar  data  sets  have  been  obtained  for  5  swirl  burner  flames,  see  Table  3.  Our  current 
focus  of  analysis  using  these  image  sets  is  the  determination  of  displacement  speed  and 
flamelet  curvature  and  characterization  of  the  reactant  flow  turbulence. 

Table  1.  V-Flames  Conditions 


Flame  # 

Mr 

uVS,° 

Un(m/s) 

Mfmm) 

1 

33.6 

0.6 

2 

0.65 

4.2 

2 

48 

1.5 

2 

0.65 

6 

3 

33.6 

1.17 

4 

0.65 

4.2 

4 

48 

2.68 

4 

0.65 

6 

5 

58.8 

0.628 

4 

0.8 

4.2 

6 

84 

1.44 

4 

0.8 

6 

7 

44.1 

0.57 

2.61 

0.7 

4.2 

8 

44.1 

1.023 

4.21 

0.7 

4.2 

9 

42 

1 

6 

0.7 

4.2 

10 

31.5 

3.16 

5.3 

0.7 

perforated  plate 

11 

48 

1.85 

3 

0.65 

6 

The  length  scale  ratio,  ////,  is  the  ratio  of  the  turbulence  integral  scale  to  the  flame 
thickness;  u  '/Si 11  is  the  turbulence  intensity  divided  by  the  unstretched  laminar  flame 
speed;  Un(m/s)  is  the  mean  flow  speed  upstream  of  the  flame  stabilizing  rod;  <f>  is  the 
equivalence  ratio;  and  finally,  M(mm)  is  the  mesh  spacing  of  the  turbulence  generating 
grid;  the  fuel  is  methane. 

Flamelet  Displacement  Speed  Measurements :  For  steady,  unstretched  premixed  laminar 
flames,  the  mass  flux  across  an  isothermal  surface  is  independent  of  the  surface  and  when 
divided  by  the  reactant  gas  density  at  the  surface  yields  the  unstretched  laminar  burning 
velocity,  Sl°.  In  a  turbulent  flame,  the  flow  is  unsteady;  the  flamelet  is  curved;  and  the 
local  mass  flux  across  isothermal  surfaces  varies  in  space  and  time.  These  variations  are 
one  measure  of  flamelet  distortion  by  turbulence.  A  related  measure  is  the  flamelet 
displacement  speed,  Sd,  that  we  define  here  as  the  component  of  the  instantaneous  gas 
velocity  relative  to  the  flamelet  and  perpendicular  to  the  instantaneous  flamelet  surface, 
as  determined  by  crossed-plane  tomography,  the  650  K.  surface. 

In  combined  CPLT  and  SPIV,  N  is  measured  along  three  parallel  measurement  lines 
(see  Figure  3)  and  Ur(x,y)  is  measured  in  a  vertical  plane.  Along  the  measurement  lines 
in  the  SPIV  plane,  N,  Uj  and  displacement  speeds  are  obtained.  To  determine 
displacement  speed,  data  on  flamelet  position,  as  marked  by  the  650  K  isotherm,  on  the 
two  measurement  lines  in  the  SPIV  plane  are  obtained  in  each  of  the  camera  images.  The 
laser  pulse  timings  of  these  images  are  staggered  and  the  flamelet  displacement  speed 
along  a  measurement  line  in  the  laboratory  reference  frame  is  obtained  by  analyzing  the 
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Results  to  Date 

Data  Taken:  Data  sets  consisting  of  CPLT  and  SPIV  images  have  been  collected  for  1 1 
different  V-flames.  The  flame  parameters  are  listed  in  Table  1,  and  the  location  of  each 
flame  on  the  Peters  regime  diagram  [9]  is  shown  in  Figure  3.  Each  data  realization 
consists  of  two  tomography  images  and  four  SPIV  images;  for  each  flame  condition 
images  from  approximately  2000  realizations  have  been  obtained  for  data  analysis. 
Similar  data  sets  have  been  obtained  for  5  swirl  burner  flames,  see  Table  3.  Our  current 
focus  of  analysis  using  these  image  sets  is  the  determination  of  displacement  speed  and 
flamelet  curvature  and  characterization  of  the  reactant  flow  turbulence. 

Table  1.  V-Flames  Conditions 


Flame  # 

l/lr 

uVS,° 

Unfm/s) 

<t> 

M(mm) 

1 

33.6 

0.6 

2 

0.65 

4.2 

2 

48 

1.5 

2 

0.65 

6 

3 

33.6 

1.17 

4 

0.65 

4.2 

4 

48 

2.68 

4 

0.65 

6 

5 

58.8 

0.628 

4 

0.8 

4.2 

6 

84 

1.44 

4 

0.8 

6 

7 

44.1 

0.57 

2.61 

0.7 

4.2 

8 

44.1 

1.023 

4.21 

0.7 

4.2 

9 

42 

1 

6 

0.7 

4.2 

10 

31.5 

3.16 

5.3 

0.7 

perforated  plate 

11 

48 

1.85 

3 

0.65 

6 

The  length  scale  ratio,  ////,  is  the  ratio  of  the  turbulence  integral  scale  to  the  flame 
thickness;  u'/Si  is  the  turbulence  intensity  divided  by  the  unstretched  laminar  flame 
speed;  Uo(m/s)  is  the  mean  flow  speed  upstream  of  the  flame  stabilizing  rod;  <f>  is  the 
equivalence  ratio;  and  finally,  M(mm)  is  the  mesh  spacing  of  the  turbulence  generating 
grid;  the  fuel  is  methane. 

Flamelet  Displacement  Speed  Measurements :  For  steady,  unstretched  premixed  laminar 
flames,  the  mass  flux  across  an  isothermal  surface  is  independent  of  the  surface  and  when 
divided  by  the  reactant  gas  density  at  the  surface  yields  the  unstretched  laminar  burning 
velocity,  Sl°.  In  a  turbulent  flame,  the  flow  is  unsteady;  the  flamelet  is  curved;  and  the 
local  mass  flux  across  isothermal  surfaces  varies  in  space  and  time.  These  variations  are 
one  measure  of  flamelet  distortion  by  turbulence.  A  related  measure  is  the  flamelet 
displacement  speed,  Sd,  that  we  define  here  as  the  component  of  the  instantaneous  gas 
velocity  relative  to  the  flamelet  and  perpendicular  to  the  instantaneous  flamelet  surface, 
as  determined  by  crossed-plane  tomography,  the  650  K.  surface. 

In  combined  CPLT  and  SPIV,  N  is  measured  along  three  parallel  measurement  lines 
(see  Figure  3)  and  Ur(x,y)  is  measured  in  a  vertical  plane.  Along  the  measurement  lines 
in  the  SPIV  plane,  N,  Uj  and  displacement  speeds  are  obtained.  To  determine 
displacement  speed,  data  on  flamelet  position,  as  marked  by  the  650  K.  isotherm,  on  the 
two  measurement  lines  in  the  SPIV  plane  are  obtained  in  each  of  the  camera  images.  The 
laser  pulse  timings  of  these  images  are  staggered  and  the  flamelet  displacement  speed 
along  a  measurement  line  in  the  laboratory  reference  frame  is  obtained  by  analyzing  the 
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stored  images  to  determine  the  flamelet  travel  distance  between  the  laser  pulses.  In  turn 
this  speed  is  used  along  with  Ur  and  N  data  to  estimate  the  normalized  flamelet 
displacement  speed  relative  to  reactants,  S^/Si0. 


-1  -0.5  0  0.5  1  1.5  2  2.5  3 

log  (I/I) 


Figure  3.  Peters  regime  diagram  showing  conditions  of  the  1 1  flames  studies.  Ka  is  the 
Karlovitz  number  based  on  //,  the  flame  thickness;  Ka§  the  Karlovitz  number  based  on  the 
high  temperature  reaction  zone  thickness,  5. 
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sd*/sL° 

Figure  4a.  Measured,  density  weighted,  flamelet  displacement  speed  relative  to  reactants. 
Data  from  the  two  measurement  lines  in  the  PIV  plane  are  shown. 
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Figure  4b.  Joint  pdf  of  curvature,  h,  scaled  by  the  laminar  flame  thickness  and  density 
weighted  displacement  speed  scaled  by  the  unstrained  laminar  flame  speed. 

The  technique  has  been  used  so  far  to  measure  Sd’/Si.0  in  a  lean,  premixed,  methane- 
air  turbulent  V-flame,  flame  1  in  Table  1.  The  displacement  speed  measured  is  that  of 
approximately  the  650  K.  isothermal  surface.  To  account  for  acceleration  within  the 
flamelet,  the  speed  is  reported  as  Sd*,  the  density  weighted  displacement  speed,  the 
product  of  the  measured  displacement  speed  and  the  density  at  650  K  divided  by  the 
reactant  density.  Flamelet  normal  data  were  also  obtained.  These  latter  data  are 
consistent  with  data  taken  in  previous  90°  CPLT  studies  of  V  and  spark-ignition  engine 
flames  [1  -  3].  The  measured  probability  density  function  of  the  flamelet  normal  is  fit  by 
a  universal  form  characterized  by  a  single  fit  parameter,  C,.  C,  is  the  same  on  the  three 
measurement  lines,  within  experimental  uncertainty. 

The  displacement  speed  results  for  this  flame  are  shown  in  Figure  4  where  it  can  be 
seen  that  the  distribution  of  Sd  /Sl°  is  peaked  with  a  mean  of  approximately  1.7±1  and 
broad  (-5<Sd/Si°<7).  The  vertical  error  bars  in  the  figure  indicate  statistical  uncertainty 
and  are  proportional  to  l/(n)' 2,  where  n  is  the  number  of  samples  in  each  histogram  bin 
used  to  construct  the  distribution.  The  horizontal  error  bars  are  estimated  by  the  rms 
error  given  in  [4,  10]  and  the  uncertainty  analysis  detailed  in  [11].  The  experimental 
uncertainty  estimates  of  the  SPIV  data  are  based  on  test  results  reported  in  [10].  Negative 
displacement  speeds  are  not  entirely  unexpected.  For  example,  they  have  been  predicted 
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to  occur  in  highly  strained  opposed  jet  flames  [12]  and  observed  experimentally  in  double 
flames  stabilized  on  an  opposed  jet  burner  [13].  In  addition  they  have  been  seen  in  a  DNS 
of  the  interaction  of  a  flame  with  a  counter-rotating  vortex  pair  [14]  and  in  a  2-D  DNS  of 
a  turbulent  methane-air  flame  [15].  The  negative  displacement  speeds  observed  in  the 
latter  calculation  are  attributed  to  regions  of  high  positive  curvature  and  extensive 
tangential  strain.  From  our  2-D  image  data  we  have  estimated  curvature  in  the  plane  of 
the  laser  illumination  sheet  and  investigated  the  correlation  between  S/  and  2-D 
curvature.  Figure  4b  shows  the  joint  PDF  of  Sd'/Si.0  and  2-D  curvature  scaled  by  the 
unstrained  laminar  flame  speed.  These  data  show  that  the  two  are  correlated  and  that 
negative  curvature  is  likely  to  increase  S/  and  the  opposite  for  positive  curvature. 

Downstream  development  of  premixed  turbulent  V-flames:  In  a  previous  study  [5]  it  was 
found  that  measures  of  flamelet  wrinkling  and  mean  turbulent  combustion  rate  increase 
linearly  with  downstream  distance  from  the  flame  holder.  A  goal  of  the  present  research 
has  been  to  use  CPLT-SP1V  data  to  explain  variations  in  the  rates  of  increase  observed 
from  flame  condition  to  flame  condition.  For  example  data  for  the  variation  of  the  PDF 
fit  parameter  versus  distance  are  shown  in  Figure  5.  In  our  work  on  this  problem  we  are 
currently  using  CPLT-SP1V  data  to  calculate  quantities  charactering  the  reactant  flow 
turbulence  and  its  downstream  development,  e.g.,  turbulent  intensity,  turbulent  integral 
and  Taylor  scales  and  components  of  the  reactant  velocities  parallel  and  perpendicular  to 
the  flamelet.  Our  goal  is  to  see  if  variations  in  these  quantities  can  explain  variations  in 
flamelet  wrinkling  and  combustion  rate  growth  between  flames. 


<;vs.  x/M 


♦  Flame  #1 
■  Flame  #2 
A  Flame  #3 
x  Flame  #4 
+  Flame  #5 

•  Flame  #6 
o  Flame  #7 
□  Flame  #9 
A  Flame  #10 
o  Flame  #1 1 


Figure  5.  Downstream  development  of  fit  parameters  for  the  eleven  flames  studied 
and  straight  line  fits  of  data. 

3-D  Curvature  determination :  The  principle  curvatures  of  a  surface  provide  a  local 
description  of  that  surface  as  the  altitude  of  the  surface  above  a  tangent  plane  of  the 
surface.  If  x,  and  y,  form  a  right-handed  coordinate  system  in  the  tangent  plane  with  its 


10 


origin  at  the  point  of  tangency  between  the  plane  and  the  surface,  the  surface  altitude 
above  the  tangent  plane  in  the  vicinity  of  the  origin  is  given  by 

hs(x,y)  =  cxx,2  +  cyy,2  ( 1 ) 

c.v  and  cy  are  the  principle  curvatures  of  the  surface  [16,  17],  and  the  quantities  that  we 
wish  to  measure.  Small  values  of  these  coefficients  mean  large  radii  of  curvature. 

In  CPLT  the  two  laser  illumination  planes  are  aligned  orthogonal  to  each  other,  and 
their  line  of  intersection,  the  measurement  line,  is  horizontal.  The  tomography  images 
can  be  processed  to  obtain  the  flamelet  curves  in  each  image  that  are  the  intersection  of 
the  650  K  isotherm  of  the  flamelet  and  the  laser  illumination  plane.  From  third  order 
polynomial  fits  of  these  curves,  flamelet  curve  tangent  vectors  are  obtained  in  each 
illumination  plane  at  the  points  where  the  flamelet  curves  intersect  the  measurement  line. 
The  cross  product  of  these  two  tangent  vectors  gives  the  flamelet  normal  vector  and 
defines  the  flamelet  tangent  plane  at  the  point  where  the  surface  intersects  the 
measurement  line. 

In  the  vicinity  of  the  measurement  line  the  flamelet  curves  in  each  illumination  plane 
can  be  expressed  in  terms  of  the  principle  curvatures  of  the  flamelet  surface.  Let  one 
laser  illumination  plane  cross  the  surface  tangent  plane  along  a  line  rotated 
counterclockwise  by  an  angle  Bp  from  the  x,  axis;  the  intersection  line  of  the  second  plane 
is  displaced  from  the  first  by  90°.  Each  illumination  plane  is  inclined  at  an  acute  angle  to 
the  tangent  plane;  <pA  for  the  first  plane  and  <ps  for  the  second.  The  specification  of  the 
angles  Bp  and  <p  for  the  two  planes  defines  their  orientation  with  respect  to  the  tangent 
plane  and  with  respect  to  x,  and  y,. 

The  flamelet  curves  in  each  illumination  plane  can  be  given  by  sets  of  points  { xp,  yp}, 
where  the  related  axes,  xp  and  yp,  form  a  right  handed  coordinate  system  in  the 
illumination  plane  with  the  jcp  axis  along  the  line  of  intersection  between  the  illumination 
and  tangent  planes.  It  should  be  noted  that  <p  can  be  determined  from  the  dot  product  of 
N  and  the  normal  vector  defining  the  orientation  in  space  of  the  laser  illumination  sheet. 

To  relate  the  measured  flamelet  curve  points  {xp,  yp)  and  the  curvatures,  one  first 
projects  these  points  onto  the  tangent  plane.  In  polar  coordinates  in  the  tangent  plane  the 
projection  relationship  is  given  by  the  following  transformation: 


(2) 


where  B,  is  the  polar  angle  of  the  projected  point  relative  to  x,  and  the  inverse  cosine  term 
is  Bm,  the  clockwise  angle  from  the  line  of  intersection  between  the  illumination  and  the 
tangent  planes  to  the  projected  point.  Substitution  of  these  points  into  Eq.  (1)  gives  hs 
where  the  surface  intersects  the  illumination  plane. 

The  altitude  of  the  illumination  plane  above  the  tangent  plane  at  points  on  the 
flamelet  curve  is  given  by 


hp  =  y„sm<p. 


(3) 


Along  the  flamelet  curve  hs  =  hp ,  and  by  combining  equations  1,  2  and  3  one  can 
obtain  a  relationship  between  cx ,  cy,  6P  and  {x p,yp}. 
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axl+flXpyp-yp+wl  =0.  (4) 

whereor  =  ( bzcx  +czcy)/ d\fi  =  2  bc(cx -cv)/a\  and  /  =  d(czcx+bzcy)/ a2. 

The  coefficients  are  a  =  tancp ,  b=cos9p,  c=sin9p  and  d=sin(p.  Here  ^  and  (p  refer  to  the  A 
or  B  illumination  plane  as  appropriate. 

CPLT  images  can  be  processed  to  give  flamelet  intersection  curves  in  the  two 
illumination  planes  as  sets  of  xp,yp,  points,  {x p,yp}.  Then  to  determine  cx, ,  cy  and  dp,  from 
these  points  first  a,  J3 and  /are  found  from  Eq.  (7)  by  minimizing  the  sum  over  {xp,yp} 
of  the  squares  of  the  left  hand  side.  Then  the  values  of  a,  ft  and  /given  above  are  used  to 
find  the  unknown  curvatures  and  angle  6P. 

To  test  this  procedure  for  finding  curvatures  a  series  of  numerical  experiments  were 
performed.  With  the  relationships  presented  above,  synthetic  or  phantom  flamelet 
intersection  curve  data  were  generated  for  arbitrary  curvature  values  and  arbitrary  values 
of  the  angles  9P,  (pu  and  (ph  defining  the  two  orthogonal  laser  illumination  planes.  These 
synthetic  data  were  analyzed  for  curvatures  and  6P. 


3-D  plot  of  plane  and  surface 


Figure  6a.  3-D  plot  of  “flamelet  surface”  and  illumination  plane  used  to  obtain  phantom 
flamelet  curve  data  shown  in  Figure  6b. 
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Figure  6b.  Phantom  flamelet  curve  data  obtained  from  intersection  of  surface  and  plane 
shown  in  Figure  6a. 
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Table  2.  Sample  Curvature  Determination  Results 


Run 

Phi(deg) 

Theta  (deg) 

cx 

cy 

# 

Plane  A 

Plane  B 

actual 

Plane  A 

Plane  B 

actual 

Plane  A 

Plane  B 

actual 

Plane  A 

Plane  B 

1 

45.00 

45.00 

45.00 

85.20 

85.20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2 

45.00 

45.00 

45.00 

45.00 

45.00 

1.00 

1.00 

3.00 

3.00 

3.00 

1.00 

3 

45.00 

45.00 

45.00 

45.00 

45.00 

1.00 

1.00 

-3.00 

-3.00 

-3.00 

1.00 

4 

18.00 

59.40 

90.00 

90.00 

0.00 

3.00 

3.00 

3.00 

1.00 

1.00 

1.00 

5 

18.00 

59.40 

90.00 

90.00 

0.00 

3.80 

3.80 

3.80 

1.20 

1.20 

1.20 

6 

18.00 

59.40 

0.00 

90.00 

0.00 

3.80 

1.20 

1.20 

1.20 

3.80 

3.80 

7 

18.00 

59.40 

120.00 

30.00 

30.00 

3.80 

1.20 

3.80 

1.20 

3.80 

1.20 

8 

36.00 

45.00 

72.00 

72.00 

72.00 

9.80 

9.80 

1.50 

1.50 

1.50 

9.80 

9 

36.00 

45.00 

178.22 

88.22 

88.22 

3.80 

1.20 

3.80 

1.20 

3.80 

1.20 

10 

36.00 

45.00 

178.82 

89.82 

89.82 

3.80 

1.20 

3.80 

1.20 

3.80 

1.20 

Sample  results  of  these  calculations  are  shown  in  Figure  5  and  presented  in  Table  2. 
For  the  calculations  two  illumination  planes  were  considered;  9pA  for  plane  A  is  specified 
while  that  for  plane  B  is  9pA  plus  ti/2.  As  noted,  different  values  of  (p  were  specific  for 
both  planes  A  and  B.  From  Table  2  it  is  seen  that  cx  and  cy  can  be  obtained  from 
curvature  data,  but  9P  is  uncertain  by  ±7t/2.  9P  specifies  the  orientation  of  the  surface  in 
rotation  about  its  normal  and  is  of  minimal  physical  significance.  Figure  5  shows  the 
intersections  of  plane  A  with  the  surface  of  row  9  in  the  table  and  the  flamelet 
intersection  curve  in  the  A  illumination  plane,  {xp,yp}.  The  dots  on  the  curves  are  the 
individual  “data”  used  for  the  fits  made  to  find  9P  and  the  principle  curvatures 

A  preliminary  test  of  the  algorithm  for  finding  3-D  curvatures  from  image  data  was 
performed  using  CPLT  images  pairs  from  flame  1  in  Table  1.  Edge  data  in  the  form  of 
third  order  polynomials  of  yp  as  a  function  of  xp  and  in  the  form  of  discrete  yp,  xp  pairs 
were  used.  The  results  to  date  have  been  inconclusive.  In  the  future  we  will  try  higher 
order  polynomial  fits  and  review  our  image  processing  algorithms  for  their  impact  on 
curvature  evaluation. 

Swirl  Burner.  In  order  to  look  at  a  different  flame  configuration  with  higher 
turbulence  levels  we  have  built  a  swirl  burner  patterned  after  one  developed  by  Robert 
Cheng  at  Lawrence  Berkeley  Laboratory.  The  flame  in  this  burner  is  free  standing,  Figure 
8  and  is  similar  to  flames  formed  behind  swirl  modules  in  gas  turbine  combustors.  Figure 
7  is  a  picture  (a)  and  a  schematic  diagram  (b)  of  the  burner,  while  Figure  8  is  a  picture  of  a 
lean  flame  stabilized  on  the  burner. 
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Figure  7a.  Picture  of  new  swirl  burner  showing  three  of  the  four  swirl  air  lines. 


Figure  7b.  Schematic  diagram  of  swirl  module  of  the  swirl  burner.  Swirl  air  is  introduced 
into  the  burner  flow  through  four  jets  as  shown. 
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Figure  8.  Short  time  exposure  picture  of  lean  («j)  =  0.7)  methane-air  flame. 

The  stable  operating  states  of  the  burner  at  different  swirl  levels,  equivalence  ratios 
and  reactants  flow  rates  have  been  determined  and  are  found  to  be  similar  to  those 
reported  by  Yegian  and  Cheng  [18]. 

Table  3.  Conditions  for  combined  CPT  -  SPIV  Measurements  in  Swril  Flames 


• 

Flame  1: 

phi  0.6 

U  =  2.6  m/s 

1 200  images 

• 

Flame  2: 

phi  0.7 

U  =  3.2  m/s 

1600  images 

• 

Flame  3: 

phi  0.7 

U  =  3.4  m/s 

1600  images 

• 

Flame  4: 

phi  0.8 

U  =  2.8  m/s 

2000  images 

• 

Flame  5: 

phi  0.8 

U  =  3.0  m/s 

2000  images 

The  data  obtained  though  these  measurements  are  on  file  but  have  not  been  analyzed 
as  of  this  date. 

Summary 

We  have  taken  sets  of  combined  CPLT-SPIV  data  for  1 1  different  methane-air  V- 
flames  covering  a  region  of  the  regime  diagram.  We  have  built  a  new  swirl  burner  and 
taken  a  similar  data  set  for  five  different  flame  conditions  from  this  burner.  From  the  V- 
flame  image  data  sets  the  turbulence  properties  of  the  reactant  flow,  flamelet  normal  data, 
displacement  speed  data  and  curvature  data  are  being  extracted.  A  method  for  finding 


displacement  speed  has  been  demonstrated  and  one  for  finding  curvature  is  under 
development.  The  swirl  burner  data  will  be  analyzed  in  a  similar  manner.  Data  of  this 
type  will  improve  our  understanding  of  how  turbulence  effects  premixed  flames  and  will 
be  of  value  for  comparison  with  the  results  of  direct  numerical  simulations  and  of  large 
eddy  simulations. 
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